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An appropriate model for the random energy landscape 
in organic glasses is a spatially correlated Gaussian field, 
generated by randomly located and oriented dipoles and 
quadrupoles. Correlation properties of energetic disorder 
directly dictates the mobility dependence on the applied 
electric field. Electrostatic disorder is significantly mod- 
ified in the vicinity of the electrode that affects injection 



properties. Correlated Gaussian field forms clusters. We 
suggest a simple method to estimate an asymptotics of 
the cluster distribution on size for deep clusters where a 
value of the field on each site is much greater than the 
rms disorder. Hopping transport in organic glasses in the 
case of high carrier density could be described in terms 
of the effective density-dependent temperature. 
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1 Introduction Most amorphous organic materials 
could be considered as organic glasses, and they belong 
to a very particular class of glasses. They are molecular 
glasses having very low concentration of intrinsic movable 
charge carriers. Hence, most organic glasses demonstrate 
measurable conductivity only after injection of carriers 
by the action of the laser pulse or strong electric field. 
At the same time, organic glasses usually have high con- 
centration of molecular dipoles and quadrupoles. Due 
to almost zero concentration of free carriers and lack of 
screening, such molecules provide long range electrostatic 
contribution to the overall energetic disorder for charge 
carriers. Usually, the resulting random energy landscape 
U(r) could be accurately approximated by the random 
Gaussian field with the typical magnitude of energetic dis- 

1/2 

order a = (U 2 ) " close to 0.1 eV Q]|2|. Long range 
sources inevitably lead to the strong spatial correlation of 
U(r): in organic polar materials (dipolar glasses) corre- 
lation function C(r) = (J7(r)Z7(0)) decays as 1/r and 
in nonpolar materials (quadrupolar glasses) it decays as 
1/r 3 [ 3 4- 1 . Long range correlations mean that sites with 
close values of random energy tend to group together and 



form clusters (cluster is defined as a set of connected sites, 
where all of them have site energy U greater than some 
boundary energy Uq) G). In this paper we consider how 
the correlated nature of organic glasses affects the impor- 
tant physical properties of these materials, namely charge 
carrier transport and injection. 

2 Charge transport in correlated Gaussian en- 
ergy landscape Correlation nature of the energy land- 
scape U(r) directly dictates major features of the charge 
transport and injection in organic glasses. For example, the 
dependence of the quasi-equilibrium mobility fx on the ap- 
plied electric field E can be understood from the following 
simple consideration (for more thorough consideration see 
Ref. 0). Suppose that the carrier is located at the bottom 
of the potential well with the energy U (0). Mobility is de- 
termined by the typical time for the carrier to reach a sad- 
dle point with the energy U(r) — eEr. That time can be 
estimated as 



t ~ CX P 



U{r) - U(0) - eEr 



kT 



(1) 
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where 1/irj is an attempt frequency, and the average time 
for the Gaussian random field U(r) is 

~U(r)-U(0) 



(t) ~ t ( exp 



to exp 




to exp 




(2) 



(7(0) - C(r) _ eEr 
~kT 



(fcT) 2 

Assuming a ID transport in the dipolar glass (DG) with 
C(r) = a 2 a/r (a is a microscopic length, comparable to 
the typical distance between neighbor molecules [ 3 , 6 1 ), we 
can calculate the critical size of the potential well that pro- 
vides the maximal escape time 



dr 



0. 



zEkT 



1/2 



(3) 



and the mobility is estimated as 



(i oc 1/ (t (r cr )) (x exp 



/ a \ 2 cr / eai? 



1/2- 



(4) 

This result provides a leading asymptotics of the exact so- 
lution of ID transport problem []6]|7]. Characteristic mobil- 
ity field dependence In fi cx v^E agrees well with the ex- 
perimental data (2j. Extensive computer simulation of the 
3D transport generally confirms Eq. (2) and only modifies 
numeric parameters in this relation. Results of the simula- 
tion may be summarizes as a phenomenological relation 
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\J eaE/a 
(5) 

where C « 0.78, and P « 2 (§). 

For the non-correlated Gaussian landscape (the Gaus- 
sian Disorder Model (GDM) (1 1), mobility is controlled by 
the carrier release from the deep states to the neighbor sites 
having higher energy, hence the shift of the carrier energy 
in the applied field lead to the linear field dependence 



In /i cx eaE/kT. 



(6) 



This estimation is in good agreement with the simulation 
data (Fig. [TJ; details of the simulation are similar to those 
described in Ref . (8) . The overall mobility dependence on 
T and E could be with reasonable accuracy described as 



1.17 - 2.05 
\kT 



eaE 



(7) 



Mobility field dependence is directly governed by the 
spatial decay of the correlation function. For example, in 
non-polar organic materials the correlation function de- 
cays as C oc 1 /r 3 and hi /j, oc E 3 / 4 EE]. In principle, 
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Figure 1 Field dependent mobility in the GDM for dif- 
ferent values of kT/a (from the top curve downward). 
Straight lines show best fits for the linear regions of the 
curves. 



the time-of-flight experiment (direct measurement of the 
carrier drift time and, hence, mobility) could be used for 
the estimation of the behavior of the correlation function 
in organic glasses. In practice, experimental complications 
(noise etc.) make very difficult (and, usually, impossible) 
to distinguish E 1 / 2 and E 3 / 4 dependences J5]]. 

3 Electrostatic energetic disorder near the elec- 
trode Structure of the organic material in the vicinity of 
the electrode is quite different from the bulk structure. 
We should expect different packing of spacious organic 
molecules, accumulation of impurities, partial degradation 
of the organic material etc. For all these reasons the ener- 
getic disorder at the interface is very different from the bulk 
disorder (typically, it is greater than the bulk disorder). 

Yet in organic glasses there is an opposite general con- 
tribution, leading to the decrease of the disorder at the elec- 
trode. Indeed, we already noted that the dominant part of 
the total energetic disorder in organic glasses is an elec- 
trostatic (dipolar or qudrupolar) disorder. The electrostatic 
energetic disorder is directly proportional to the disorder 
in the spatial distribution of electrostatic potential, gener- 
ated by molecular dipoles or quadrupoles. In organic layers 
sandwiched between conducting electrodes this spatial dis- 
tribution must obey a boundary condition at the electrode 
surface: here the potential should be a constant. Thus, at 
the electrode there is no electrostatic disorder at all, irre- 
spectively to how disordered is the material in the bulk. 
This means that the magnitude of the dipolar or quadmpo- 
lar disorder increases while going away from the interface, 
asymptotically reaching its bulk value. 
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Some decrease of the disorder at the surface of organic 
material is inevitable for any model of organic glass (just 
because there are more neighbor molecules in the bulk of 
the material), but the magnitude of the effect in the case 
of electrostatic disorder is much greater than in the case of 
short range interactions. For example, for the simple model 
of the interaction with the nearest neighbors only and sim- 
ple cubic lattice we have of urfacc = 5/6cr^ ulk , while for 
the dipolar disorder CTg Urfaco ~ 0.3cr^ ulk and rms disorder 
depends on the distance z from the interface as 



'(*) 



' bulk 



1-^(1- 

2z 



-2z/a 



ao = Aa, 



(8) 

and for a simple cubic lattice A = 0.76 flOl . Spatial corre- 
lations at the interface do differ too; this phenomenon has 
no analogue for the short range disorder. A direct calcu- 
lation of the correlation function C(r) near the interface 
gives 



C(zi,z 2 , p) = cr bu i k a 



1 



1 



(9) 



where r ± = p 2 + (z\ ± z^) 2 and p is a 2D vector oriented 
along the interface plane ifTTI . Hence, at the interface the 
dipolar glass is much less correlated in comparison with 
the bulk: C(zi, z%,p) oc Z1Z2I P 3 for P ^ z%, and clus- 
ters are elongated perpendicular to the interface. 

Decrease of the disorder at the interface and change 
of the spatial behavior of the correlation function have 
very important implications for the charge injection. In 
the absence of these phenomena, injection current in or- 
ganic glasses demonstrates formation of channels where 
the current density is much greater than the average den- 
sity lfl2l . Such channels originate from particular spots at 
the interface, where clusters of sites with low energy facil- 
itates injection. Reduction of the disorder at the interface 
and modification of the spatial behavior of the correlation 
function lead to the more uniform distribution of the in- 
jection current over the electrode and dramatically reduce 
current channeling. This decreases local overheating in a 
device and improves its performance. 

4 Distribution of cluster size for deep clusters A 

common feature of any random medium is the formation 
of clusters. One of the basic characteristics of the cluster 
statistics is the cluster distribution on size (or cluster num- 
bers) n s , here s is the number of sites in a cluster. 

Analytical results for n s are scarce, even for the sim- 
plest case of non-correlated disorder. The reason for the 
scarcity is obvious: it is difficult to take into account var- 
ious shapes of clusters. Most results in this area were ob- 
tained using scaling arguments with subsequent testing of 
their validity with extensive computer simulation II13II14I 

El 

It turns out that a very simple calculation of the distri- 
bution Py (Uq) of the average value of the correlated Gaus- 
sian field Uq in finite domain having volume V provides a 



very accurate estimation of the leading asymptotic for the 
distribution of large clusters with s ^> 1 and Uq ^ a. 

The distribution Py (Uq) for the Gaussian random field 
can be calculated exactly 1 16| 



Pv(U ) = 



V 



V2ttK 



= ex P 



u 2 v 2 

2K 



(10) 



K= ( drdriC(r-r-i). 
Jy 



For a spherical domain with radius Rq in the dipolar glass 



K = ^-Aaa 2 R\ oc V b '\ 
15 



(11) 



while for the non-correlated field K oc Rq oc V lfl6l . 

We may expect that Eqs. (11011 U give a reasonable esti- 
mation for the number n s of the true clusters, i.e. domains, 
where U(r) > Uo everywhere (assuming V = a 3 s), at 
least for the leading term of the asymptotic dependence of 
n s on s (the very use of the continuous model of the ran- 
dom medium suggests that our consideration is valid only 
for s S> 1). If so, then in the dipolar glass 



Inn, oc 



B 



4A(36tt) 1 /3 



0.34. 



(12) 

This hypothesis agrees well with the simulation data Ifl6l . 
For the non-correlated Gaussian field it gives a well-known 
exact asymptotics lnn s oc —U 2 s/2cr 2 IfTTI . 

5 Charge transport of interacting carriers in 
Gaussian glasses If charge carrier density is not very 
low, we cannot neglect the Coulomb interaction between 
carriers. The relevant interaction strength parameter is 
A = e 2 n 1 l' i I ea, where e ~ 2.5 — 3 is a dielectric constant 
of the organic material and n is a carrier concentration. 
The maximal value of A for n = 1/a 3 in organic glasses is 
Amax — 5. This interaction provides an additional compli- 
cation to carrier dynamics. Yet spatial correlation manifests 
itself even in the case of high carrier density. 

Indeed, a general tendency for the DG model is that 
transformation of the mobility curve with the increase of 
average fraction of the occupied sites p = na 3 resembles 
the corresponding transformation of the curve with the in- 
crease of T (compare insets in Fig. [2] for example, with 
Fig. 1 in Ref. (8J): with the increase of p mobility becomes 
greater and the slope of p,(E) curve becomes smaller. This 
is not the case for the GDM: here only mobility curve 
moves upward but the slope remains approximately con- 
stant. Details of the simulation can be found in Ref. Ifl8l 

This difference could be easily understood. It was 
noted that the field dependence of [i in the GDM is gov- 
erned by the carrier escape from deep states to the nearest 
sites having much higher energy, and the field-induced 
shift of site energies leads to Eq. ©. Random charge 
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Figure 2 Dependence of the effective temperature on the 
carrier density p for kT/a = 0.3 (■) and kT/a = 0.4 
(□) for the DG model (top plot) and GDM (bottom plot), 
correspondingly. Insets show the curves n{E) for kT/a = 
0.3 and different p: 0.0016, 0.008, 0.016, and 0.032, from 
the bottom curve upward, correspondingly. For the DG T c ff 
was calculated by fitting the mobility data to Eq. ©, while 
for the GDM the mobility data was fitted to Eq. (01. Lines 
are provided as guides for an eye. 



distribution provides a smooth random energy landscape 
superimposed on the intrinsic disorder, but typical addi- 
tional variation of energy at the scale a is negligible for 
small p. Hence, estimation (O remains valid and the slope 
of the mobility curve does not depend on p. 

Situation in the DG model is different: here mobility 
field dependence is governed by the carrier escape from 
critical clusters, as described in Section If we increase 



the density of carriers, then at first they fill these critical 
traps, because the release time is maximal here. Hence, 
transport of more mobile carriers is governed by clusters 
with the size that differs from r cr (it is smaller). This means 
that the effective critical size r„ depends on p. According 
to Eq. (0, this is equivalent to the dependence of the effec- 
tive temperature T e g on p: T e g should grow with p. This 
conclusion is in good agreement with Fig. [2] while for the 
GDM T c ff does not depend on p and is very close to T, for 
the DG model T e g monotonously grows with p. 

6 Conclusion In organic glasses the dominant part 
of the random energy landscape is formed by the elec- 
trostatic contributions from randomly located and oriented 
dipoles and quadrupoles. Long range electrostatic contri- 
butions automatically provide long range spatial correla- 
tions of random energies. These correlations directly dic- 
tate the form of the mobility field dependence for the hop- 
ping charge transport in amorphous organic materials, even 
in the case of high carrier density, where carrier-carrier in- 
teractions cannot be neglected. Electrostatic disorder is re- 
duced at the interface with the electrode, thus significantly 
modifying the injection properties. A simple but effective 
procedure for the calculation of cluster distribution on size 
for a correlated random Gaussian field is suggested. 
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